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A class of solutions of the gravitational field equations describing vacuum spacetimes 
outside rotating cylindrical sources is presented. A subclass of these solutions corresponds 
to the exterior gravitational fields of rotating cylindrical systems that emit gravitational 
radiation. The properties of these rotating gravitational wave spacetimes are investigated. 
In particular, we discuss the energy density of these waves using the gravitational stress- 
energy tensor. 
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1. INTRODUCTION 



The main purpose of this paper is to present a class of cylindrically symmetric vacuum 
solutions of the gravitational field equations representing rotating gravitational waves and 
to study some of the physical properties of such waves. The spacetimes under consideration 
here are not asymptotically flat in general; therefore, the concepts of energy, momentum, 
and stress do not make sense in the standard interpretation of general relativity. Never- 
theless, it is possible to introduce a local gravitoelectromagnetic stress-energy tensor via a 
certain averaging procedure [1]. This gravitational stress-energy tensor provides a natural 
physical interpretation of the Bel-Debever-Robinson tensor that has been used frequently 
in numerical relativity [2]. In this work, we study the local gravitational stress-energy 
tensor for free rotating gravitational waves. 

Rotating cylindrically-symmetric gravitational waves were first discussed in 1990 [3, 4]. 
The investigation of these solutions was motivated by Ardavan's discovery of the speed- 
of-light catastrophe [5] and its implications concerning gravitation [6]. Previous work is 
generalized in the present paper and an extended class of rotating gravitational wave space- 
times is analyzed. We find that the solution investigated previously [1, 2] is special, since 
it is the only member of the extended class studied here that represents the propagation 
of free rotating gravitational waves. 

In this paper, we consider Ricci-flat spacetimes that are characterized by a metric of 
the form 

-ds 2 = e 2 ^(-dt 2 + dp 2 ) + p 2 e~ 2 ^(udt + #) 2 + e 2 ^dz 2 , (1) 

in cylindrical coordinates (p, </>, z). Here 7,//,^ an d to are functions of t and p only; 
moreover, the speed of light in vacuum is set equal to unity except where indicated other- 
wise. The spacetime represented by equation (1) admits two commuting spacelike Killing 
vectors d z and d^. Though d z is hypersurface-orthogonal, is not and this fact implies 
that the isometry group is not orthogonally transitive. Instead of the variables t and p, it 
is convenient to express the gravitational field equations in terms of retarded and advanced 
times u = t — p and v = t + p, respectively. The field equations then take the form 
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(fJ>1pv)u + {p4u)v = 0, (2) 

f*uv - ^" 3e27 = 0» ( 3 ) 

u v - u u = l(x~ 3 e 21 , (4) 

lu = 7) (Vuu + Z^l), (5) 

Iv = ^— (fivv + 2//V«), ( 6 ) 

Lfl v 



where tjj u = dip/du, etc. Here / is a constant length characteristic of the rotation of the 
system. Using equation (3), it is possible to eliminate 7 from equations (5) and (6); then, 
we obtain the following equations 

Wu = 3 /"u - A»A»u« + WuHuuviHuv)' 1 (7) 

2// 2 ^ = 3^ - + HUvHuwiHuvY 1 - (8) 

The integrability condition for this system — i.e. = — results in a nonlinear 

fourth order partial differential equation for \x. Alternatively, one could obtain the same 
equation for fx by combining equations (2), (7), and (8), which shows the consistency of 
the field equations (2) — (6). It is important to notice that fi cannot be a function of 
u or v alone, since this possibility would be inconsistent with equation (3). The partial 
differential equation of fourth order for fx is, however, identically satisfied if \x is a separable 
function, i.e. fx = a(u)f3(v); this leads, in fact, to the rotating waves discussed earlier [3, 4]. 
Here we wish to study the general solution of the field equations. 

In section 2, we discuss the general solution of the field equations. We find that there 
are two possible classes of solutions: The first class corresponds to the stationary exterior 
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field of a rotating cylindrical source, while the second class appears to represent a mixed 
situation involving rotating gravitational waves. In fact, a subclass of the latter solutions 
describes the exterior fields of certain rotating sources that emit gravitational radiation; 
these solutions approach the special rotating wave solution [3, 4] far from their sources. 
Indeed, the only solution of the second class corresponding to a pure gravitational wave 
spacetime is the special solution [3, 4] that is further discussed in section 3. In section 4, 
some aspects of the energy and momentum of the special rotating gravitational waves are 
discussed using the gravitoelectromagnetic stress-energy tensor developed in a recent work 
[1]. The appendices contain some of the detailed calculations. 



2. SOLUTION OF THE FIELD EQUATIONS 

To solve the field equations (2) — (6), let us introduce the functions U, V, and W by 

U = Hu^u ~^uu, (9) 

V = [l v ^ v - -fJL vv , (10) 

3 , . 

W = W uv + -fJL uv , (11) 

and rewrite equations (5) and (6) as 

^ u = e(fiU) 1 / 2 , ^ v = e{^V) 1/2 , (12) 

where the symbols e and e represent either +1 or —1 (i.e., e 2 = e 2 = 1). Using equation 
(3), it is straightforward to show that 



HU V = Lt u W, fj,V u = fjLyW. (13) 

Let us now combine relations (12) and (13) in order to satisfy equation (2); the result is 
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e U-^ifivU + fi u W) + e V- 1 ' 2 ^ + fi v W) = 0. 
This equation can be written as 



(14) 



Uiji u V + ii v Wf = V{fi v U + VuW) 2 , (15) 

which holds if either W 2 = UV or /x£ V = /j, 2 U; that is, equation (15) can be factorized 
into a fourth order equation and a third order equation, respectively. 

Let us first consider the fourth order equation W 2 = UV. It follows from this relation 
and equation (12) that W = ±fj,ip u ip v . This expression for W along with U = flip 2 and 
V = fiifi 2 can be inserted into equation (13) to obtain relations that when combined with 
equation (2) in the form 

tyipuv = -(Vuipv + Vvipu) (16) 

imply 

Mv = o. (17) 

If ip is a function of either u or v alone, then equation (2) implies that ip must be a constant. 
The spacetime represented by a constant ip turns out to be a common special case of the 
solutions discussed in this section and is flat. 

Let us next consider y? u V = /j, 2 U . It follows from equation (3) that 

^ 3 /i M ^ 1 [i uvu (18) 

2 fj> 2 fi U V 

2 fj> 2 fj> U v 

which together with equations (9) and (10) imply that y? u V = ji 2 U is essentially equivalent 
to 



2 ~~ 2 ' V l 

H>v A*d A*u A*m 
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Equation (20) can be written as 



(\n/i v ) uv = (\n/i u ) vu . (21) 

Let f(u) and g{v) be arbitrary functions of their arguments and consider the transformation 
(n, v) — > (x, y) such that 

x = f(u) + g(v), (22) 



y = f(u)-g(v), (23) 

f u 7^ 0, and 7^ 0; then, the general solution of equation (21) is that \x should not depend 
on y, i.e. 

H = n(x). (24) 

Let us now proceed to the calculation of ip. To this end, let us consider a function 
C(x) given by 

£ 2 (z) = i(V 2 -^'' + ^^), (25) 

where a prime indicates differentiation with respect to x. Equations (7) and (8) can now 
be written as 



^ 2 u = f 2 C 2 , fy 2 v =g 2 C 2 . (26) 



Furthermore, 
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Combining equations (26) — (28), we find that either ip is purely a function of x given by 
fj, 2 (dip / dx) 2 = C 2 (x) or ip is purely a function of y given by /j, 2 (dip / dy) 2 = C 2 (x). These 
cases will now be discussed in turn. 



Case (i): ifj = ip(x) 

It follows from equation (2) that in this case 

~r~ ~ t~t i (29) 
dx n(x) 

where C is an integration constant. Thus /j(x) is determined by C 2 (x) = C 2 . To solve this 
differential equation, let X = p! and note that 

-2 _ 1 V 2 ( 2 dx \~ X d ( 3 dX \ 



C -2 X \^) m\fni ■ (30) 



Moreover, let S = ji dX/dfi; then, C (x) = C can be written as 



ds _ 2 fc 

dX~ \X 



= 2 ( — 2 ~ 1 ) , (31) 



which can be integrated to give S = —2(D+X + C 2 /X). Here D is an integration constant. 
It follows from this result that fi(x) can be found implicitly from dfi = X(fx)dx, where 

^(-/ jTO+c* )- < 32 > 

We will show that for C = 0, ip is constant and this solution turns out to be flat. If C ^ 0, 
however, we have a new solution. 

The spacetime metric in this case can be written in a form that depends only on x. 
That is, it is possible to show — by a transformation of the metric to normal form - 
that the general solution in Case (i) has an extra timelike or spacelike Killing vector field 
d y . To this end, let us consider the coordinate transformation (t,p,<f>,z) — > (T, R, <E>, Z), 
where f(u) = (T + R)/2, g(v) = (-T + R)/2, = $ + H(T,R), and z = Z. In this 
transformation H(T, R) is a solution of the partial differential equation 

_ ( , + G) §* +{F - G fJL + l i(F + G) ^ , ,33) 
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where F = f u 1 , G = g v 1 7 R = x, and T = y. To prove our assertion, we will show that 
under this coordinate transformation the spacetime metric in Case {%) takes the form 



-ds 2 = P(R)(-dT 2 + dR 2 ) + Li 2 (R)e- mR) MR)dT + d<S>} 2 + e 2 ^ R) dZ 2 , (34) 

which is clearly invariant under a translation in T thus implying the existence of a Killing 
vector field <9t- If P(R) > 0, then this metric could represent the stationary exterior field 
of a rotating cylindrical configuration. It follows from a comparison of the metric forms 
(1) and (34) that 

p(fl) = 4" 3 % e " 2 * < ' !, ' (36) 



r) FT 1 

w = - 1 oj(F-G) + Q(R), (36) 



r) M 1 

- = - r(F + G) . ,37) 

Equations (36) and (37) can be used in equation (33) to show that Q(R) = 2l~ 1 dfi/dR. It 
remains to show that the integrability condition for equations (36) and (37), i.e. 
d 2 H I dRdT = d 2 H/dTdR, is satisfied. It turns out that this relation is indeed true, since 
it is equivalent to the field equation (4) for ui. 

Finally, let us consider the case C = or, equivalently, ip = constant. It follows from 
equation (31) that in Case (z) 



d 2 \i 2 



2 

dR +D dR +C 



(38) 



dR 2 /I 

and one can show explicitly that the metric form (34) is flat once C = 0. For C ^ 0, the 
general solution of Case {%) is not flat; however, its physical properties will not be further 
discussed in this work, which is devoted to rotating gravitational waves. 



Case (ii): if) = ip(y) 
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It follows from (dip / dy) 2 = C 2 (x)//j 2 that ip must be a linear function of y, since the 
left hand side of this equation is purely a function of y while the right hand side is purely 
a function of x; therefore, each side must be constant. Thus ip = ay + b, where a and 
b are constants, and C 2 = a 2 /i 2 . It turns out that equation (3) is identically satisfied in 
this case. Using equation (30) , the differential equation for \i can be expressed in terms of 
X = \i' as 

x2 d ( 3 d_X\ =2(? ^dX 
dfx \ d[i ) d\x 

If a = 0, ip is constant and we recover the same flat spacetime solution as in Case (i) 
with C = 0. Therefore, let a ^ and consider a new "radial" coordinate r given by 
r = (a 2 ^ 2 ) -1 ; then, 

r 2 X 2 — r + — = 0. 40 
dr l dr 

The solution X = constant is unacceptable, since it implies that 7 = —00. However, there 
is another exact solution that is given by 

3 ^ 1/2 



X = ± ( - rj , (41) 

which turns out to correspond to the special rotating gravitational waves [3, 4] that are 
the subject of the next section. It is possible to transform equation (40) to an autonomous 
form; to this end, let us define A and B such that A = 2/(3rX 2 ) and B = -4X r /(3X 3 ), 
where X r = dX/dr. It can then be shown that equation (40) is equivalent to 

d_B _3 B(B-A 2 ) 

dA 2 A(B - A) ' 1 ' 

where (A, B) = (1,1) represents the special solution (41). This special solution is an 
isolated singularity of the nonlinear autonomous system (42) and the behavior of charac- 
teristics near this point indicates that (1,1) is a saddle point. 

Let us note here certain general features of X(r), which is the solution of equation 
(40). If X(r) is a solution, then so is —X(r). Moreover, equation (40) can be written as 
(X 2 ) r = —2X 2 /(r 2 X 2 ), which indicates that for X 7^ constant the absolute magnitude of 
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the slope of X(r) monotonically decreases as r increases. If X(r) has a zero at r ^ 0, 
then the behavior of X(r) near r$ is given by 



2 r x 



X(r) = ±J— 1 

r V^o 



1 f r \ 3 / r 



(43) 



for r > r$. These results are illustrated in figure 1. 

Once a solution X(r) of equation (40) is given, one can find a solution of the field 
equations in Case (ii) . It turns out that in general such a solution is not a pure gravitational 
wave. The only exception is the special solution (41). To demonstrate this, let us first 
consider a transformation of the spacetime metric (1) to the normal form appropriate for 
Case (ii). The first step in this transformation will be exactly the same as that given 
in Case (i) and equation (33), i.e. (£, p, 0, z) — > (T, i?, $,Z). Next, we change the radial 
coordinate to r and let (T, R, Z) — > (i, r, Z), where i = 2aT and Z = la 2 exp(2b) Z . 
In terms of the new system of coordinates, the spacetime metric may be expressed — up 
to an overall constant factor — as 

- ds z = e -i (^-X 2 dP + ^P)+ ^-(iXdi + d$) 2 + e { dZ 2 , (44) 

where I = (la) -1 is a constant and X(r) is a solution of equation (40). There are four 
curvature invariants for this Ricci-flat spacetime and it is possible to show that two of these 
are identically zero. The other two also vanish for the special solution (41); therefore, this 
special case corresponds to free gravitational waves [7] . Specifically, let 

h = R^upaR^ ~ iR^paR*^, (45) 

h = R^R^RJ" + iR^R^R^", (46) 

be the complex invariants of spacetime. We find that for the metric form (44) these are 
real and are given by 

7i = -e 2i (rX 4 X 2 )-\l - 3rX 2 - 2r 2 XX r - r 3 X 3 X r + r 4 X 2 X 2 ), (47) 
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I 2 = --e 3i (rX 5 X^)- 1 (l - 2rX 2 - 2r 2 XX r + 2r 3 X 3 X r + r 4 X 2 X 2 ). (48) 

It is important to note that Ii and I 2 vanish for the special solution X = ±(3r/2) -1 / 2 , 
indicating that it describes the propagation of free gravitational waves. Moreover, all other 
solutions — which do not describe free gravitational waves — are singular for t — > oo, since 
7i and I2 both diverge in the infinite future. Thus all such solutions evolve to states that 
are ultimately singular. The singular nature of the special solution (41) has been discussed 
previously [3,4]. 

It is clear from equation (40) and figure 1 that as r — > 0, X(r) — > ±(3r/2) -1 / 2 for 
a subclass of the spacetimes under consideration here. Such solutions correspond to the 
exterior field of a rotating and radiating cylindrical source such that very far from the 
symmetry axis the metric nearly describes free rotating gravitational waves given by the 
unique special solution (41). Indeed, inspection of the metric form (44) reveals that the 
circumference of a spacelike circle about the symmetry axis and normal to it is proportional 
to r -1 / 2 at a given time i; therefore, as r — > the metric form (44) describes the asymptotic 
region very far from the axis. Moreover, it is simple to check that I\ — > and I2 — > as 
r — > for the solution given in figure 1 with X r (r = 1) = —1. 

Finally, let us note that the unique special solution corresponding to free gravitational 
waves (41) can be written in terms of \i and x as dji / dx = ±(2/3) 1 / 2 a / u, which can be easily 
solved to show that fx depends exponentially upon x = f(u) + g(v). This means simply 
that fi can be written as \i = a(u)f3(v), where a and (3 are arbitrary functions. The next 
section is devoted to a discussion of the properties of this solution beyond what is already 
known from previous studies [3,4]. 

3. FREE ROTATING GRAVITATIONAL WAVES 

Free nonrotating cylindrical gravitational wave solutions of Eintein's equations were 
first discussed by Beck [8]. These solutions have been the subject of many subsequent 
investigations [7,9,10]. The solutions considered in this section can be interpreted in 
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terms of simple cylindrical waves that rotate [4] . 

The special solution (41) for free rotating gravitational waves is given by 



li = a(u)p(v), (49) 
*=\j\^ (50) 

(51) 

u v - u> u = ja u (3 v . (52) 

It turns out that in this case equation (2) is equivalent to the scalar wave equation for tfj 
in the background geometry given by equation (1); therefore, the function ip — which is a 
mixture of ingoing and outgoing waves according to equation (50) — has the interpretation 
of the scalar potential for the free rotating gravitational waves. The solution (49) — (52) 
cannot be thought of as a collision between outgoing and ingoing gravitational waves, since 
the field equations do not admit solutions for which p = p(u) or p = p(y). That is, there 
is no purely outgoing solution just as there is no purely ingoing solution. 

The spacetime given by equations (49) — (52) is singular. In fact, the analysis of 
the corresponding curvature indicates that moving singular cylinders appear whenever 
a,/3,a u , or j3 v vanishes. It is interesting to consider the nature of the symmetry axis for 
rotating waves, since in these solutions the axis does not in general satisfy the condition 
of elementary flatness. If for an infinitesimal spacelike circle around the axis of symmetry 
the ratio of circumference to radius goes to 2n as the radius goes to zero, the condition of 
elementary flatness is satisfied for the axis under consideration. In our case, this means 
that p 2 /{e 2l p 2 ) — > 1 as p — > 0. For simple cylindrical waves (i.e., Beck's solution) we 
have p = p and hence the condition of elementary flatness is that 7 — > as p — > 0. In 
general, Beck's fields can be divided into two classes: the Einstein-Rosen waves [9] and 
the Bonnor- Weber- Wheeler waves [10]. In the former class, the axis does not satisfy the 
condition of elementary flatness; in fact, the axis is not regular either. It is a singularity 
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of spacetime and is therefore interpreted as the source of the cylindrical waves which are 
otherwise free of singularities. In the latter case, the axis does satisfy the condition of 
elementary flatness and is, moreover, regular. The waves presumably originate at infinity: 
Incoming waves implode on the axis and then move out to infinity with no singularities 
in the finite regions of spacetime. It is therefore clear that no caustic cylinders appear 
regardless of the nature of the axis. They do appear, however, when the waves rotate. We 
are therefore led to regard the appearance of moving singular cylinders in our solution as 
being due to the rotation of the waves. It is important to emphasize the absence of a direct 
causal connection between the violation of elementary flatness at the axis and the presence 
of singular cylinders: the axis is static while the singularity is in motion, and the condition 
of elementary flatness involves the gravitational potentials {g^ v ) while the singularity of 
the field has to do with the spacetime curvature. In fact, it is possible to find instances of 
exact solutions for which the axis is elementary flat but singular [11]. 

In spacetime regions between the extrema of a and /?, the solution (49) — (52) can 
be reduced to a normal form in which a and (3 are linear functions of their arguments [3] . 
With further elementary coordinate transformations, the normal form can be reduced to 
two special solutions which are of interest: (a) a = cm, (3 = av and (b) a = au, j3 = —av. 
Here a = ±(/o) -1 ^ 2 , where Iq is a constant length whose introduction is necessary on 
dimensional grounds. The metrics for these cases are, respectively, 



-* 2 = w {uvf + * 2) + M [uvf + # ) 

and 



III 2 



u 



ds z = --^{-uvf [--) (-dt 2 + dp 2 ) + ^(uv 



f 2 
'-o 



u 



T-jPdt + 



+ 



(-1) 



V6 



dz 2 . 



(54) 
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To ensure that the spacetime metric is real, case (a) must be limited to (u > 0, v > 0) or 
(u < 0,v < 0). Similarly, case (b) must be limited to (u > 0, v < 0) or (u < 0, v > 0). 
In either case the hypersurfaces u = and v = are curvature singularities. Since the 
laws of physics break down very close to these surfaces, it appears that the consideration 
of boundary conditions across such surfaces would be without physical significance. In 
case (6), p is a temporal coordinate and t is a spatial coordinate. The transformation 
(t — > p, p — > t), or equivalently (it — > —it, u — > u), brings the metric to the form 

+ (^)^& 2 , (55) 

which reduces to case (a) with a further coordinate transformation — > + 8tp/(lal). 
Therefore, only the normal form in case (a) will be considered in the rest of this paper. 

The geodesic equation for free rotating gravitational wave spacetimes is discussed in 
appendix A. 



4. THE GRAVITATIONAL STRESS-ENERGY TENSOR 

In a recent paper [1] , a local gravitoelectromagnetic stress-energy tensor T M „ has been 
defined for Ricci-flat spacetimes via a certain averaging procedure in a Fermi frame along 
the path of a geodesic observer. That is, 

L 2 

M " ~ 12-k ^ vpa (o) (o) 
L 2 

= T^ T ^(°)(°) ' ( 56 ) 

where L is a length-scale characteristic of the field under consideration, T^ upa is the Bel- 
Debever-Robinson tensor [7, 12] defined by 

T^upa = 7^ (Rfj£pcRv£ + p) ~ Y^9p,u9paRai3 1 sR at3 ' lS ', (57) 
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and = dx^/dr is the vector tangent to the timelike path of the observer. The grav- 
itational stress-energy tensor is expected to provide approximate measures of the average 
gravitoelectromagnetic energy, momentum, and stress in the neighborhood of the observer. 
There does not seem to be any direct connection between T^ v and the Landau-Lifshitz 
pseudotensor; this point is further discussed in appendix B. 

It would be interesting to study the energy of rotating gravitational waves. However, 
the spacetime metric for the rotating waves (44) is given explicitly only for the case of the 
exact solution (41) of the differential equation (40). This unique solution is approached 
asymptotically by the rotating wave solutions very far from the symmetry axis. For the 
sake of simplicity, we therefore restrict attention to the special solution (41) that, according 
to the results presented in section 3, corresponds to the free rotating waves. 

To determine the energy density of this radiation field measured by a geodesic observer, 
it is first necessary to consider a solution of the geodesic equation discussed in appendix 
A. Let us therefore choose a "radial" geodesic of the free rotating waves given by 

t= l j(C 1/p + C 1/q ), (58) 




z = z , (61) 

where we have used equations (All) — (A14) of appendix A together with the metric in 
normal form (53), i.e. a(u) = u/l^ 2 and (3(v) = v/l^ 2 , and we have fixed the constants 
in these equations such that t\ = r 2 = lo and 

2 C1 = ll~ p/ \ 2c 2 = /J"* 72 . 

Note that the scaling parameter Iq is therefore related to the intrinsic rotation parameter 
I via 2/o = V^l. Moreover, </> and zq are constants and £ is defined by 
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C = l-r// . (62) 

The observer under consideration here starts at the symmetry axis p = 0, which is regular 
at t = 0, and returns to it at r = Z when it is singular. The trajectory of this observer 
p{4>) is depicted in figure 2. 

Imagine now an orthonormal tetrad frame A^ s that is parallel propagated along this 
timelike worldline. It is simple to work out explicitly two axes of the tetrad: the time axis 
and the spatial axis parallel to the symmetry axis. These are given by A^ = dx^/dr using 
equations (58) — (61) and 

A (3) = C 3/5 C (63) 

respectively. It follows from the projection of T^ vpa given by equation (57) on these axes 
that along the geodesic the gravitational radiation energy density is given by 

36 L 2 

T(0 ^ 0) = 6257T (l Q -rr ( ] 

while the energy flux along the symmetry axis vanishes, ^o)(3) = 0, as expected. There is, 
however, radiation pressure along the z-axis and is given by 

3 I? 

T(3)(3) = 1257 Tk^W (65) 
The energy density and the pressure measured by the geodesic observer both diverge at 
the singularity r = / . Note that T( 3 )( 3 )/T( )(o) = 5/12, which is less than unity as would 
be expected for the ratio of pressure to density. The characteristic length-scale associated 
with rotating gravitational waves is /; therefore, the constant length L in expressions (64) 
and (65) could be chosen to be simply proportional to / = 2/ /\/5- 



APPENDIX A: GEODESICS IN ROTATING WAVE SPACETIMES 



Starting with the metric form given by equation (1), let 
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2 

C g = -\ = \[e^{-i 2 + P 2 ) + A-^M + </>) 2 + e^z% (Al) 

where i = dt/dX, etc. Since the Lagrangian (Al) does not depend upon 4> and z, there are 
two constants of the motion p^ and p z given by 

P^ = ^f=p 2 e-^(ui + <P) 7 (A2) 

Moreover, we have 
and 

Let us now focus attention on "radial" geodesies such that the momenta associated 
with azimuthal and vertical motions vanish. It is then convenient to write the equations 
of motion in terms of radiation coordinates u and v. We find that equations (A2) — (A5) 
reduce to 

and 

1 Ms) ~ n f x % <-) 



which imply that 



e 27-2?/> du dv 



dX dX 

is a constant along the path. This constant vanishes for a null geodesic, so that "radial" 
null geodesies correspond to u = t — p = constant or v = t + p = constant, where A is the 
affine parameter along the path. 
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We are mainly interested in timelike "radial" geodesies, hence we set A = r, where 
r is the proper time along the path. Then, with U = du/dr and V = dv/dr, we have 
UV = e 2 ^" 27 . It follows from (A6) and (A7) that U V V 2 = V U U 2 . The general solution is 
given by 

du 1 f dS\ 1 , .„ s 

rfw ^v 1 , { A9) 



dr 2\dv J 
where S(u,v) is any solution of the differential equation 

* ££ = 

This is, in fact, the Hamilton- Jacobi equation for "radial" timelike geodesies, and S along 
the path differs from the proper time by a constant. 

It is simple to illustrate a class of solutions of the eikonal equation (A10) via separation 
of variables, i.e. 

S(u, v) = - Cl a p (u) - c 2 /3 q (v), (All) 
where p = 4 — \/6, q = 4 + and c\ and c 2 are constants such that 

Cl c 2 = ± (A12) 
Equations (A8) and (A9) now have solutions 

a(u) = (^r) 1/P - (413) 
«°> = ( ! ^ I ) 1/ '' iAU) 

where T\ and r 2 are constants of integration. More explicitly, let us consider the normal 
form of the metric with a(u) = l l ^ 2 u and j3(v) = l 1 ^ 2 v as employed in section 4. The 
geodesic would then hit the null singular hypersurface u = (or v = 0) at r = t\ (or 
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t = r 2 ). The singular nature of these moving cylinders can be seen from the fact that the 
geodesic can not be continued past them since then u or v would become complex. 



APPENDIX B: THE LANDAU-LIFSHITZ PSEUDOTENSOR IN RIEMANN 
NORMAL COORDINATES 

The energy of a gravitational field — if it can be defined at all — is nonlocal ac- 
cording to general relativity. On the other hand, the Bel-Debever-Robinson tensor is 
locally defined. A connection could perhaps be established between these concepts if the 
energy-momentum pseudotensor of the gravitational field is expressed in Riemann normal 
coordinates about a typical event in spacetime. 

Let x M be the Riemann normal coordinates in the neighborhood of some point ("ori- 
gin") in spacetime; then, 

9ixv = t}ixv ~ ^Rp.au/3 X a X P H , (Bl) 

K P = -\(R%* + R%v«)x° + ■■■■ m 

The Landau-Lifshitz pseudotensor is quadratic in the connection coefficients by construc- 
tion; therefore, t^~ L is — at the lowest order — quadratic in Riemann normal coordinates. 
Hence, 

c 4 

^"-^ = 144^G eM " Q/3 + '"' (S3) 
where Mi , Q/ 3 is symmetric in its first and second pairs of indices by construction and is 
given by 



1 7 

= ~2^R P ^ a ^ UCJ PP ^ ^pRwpa) "I" ~^(Rlip<ja.Rv p Rp-pcrpRJ' a ) 

-lv^v^R P ^sR prTK5 . (B4) 
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This expression should be compared and contrasted with equation (57) that expresses the 
Bel-Debever-Robinson tensor in a similar form. There is no simple relationship between 
®iiva[3 and T Mi , Q/ 3; however, one can show that 

in Riemann normal coordinates. Here K is the Kretschmann scalar, i.e. K = R^ upa R^ vpa . 
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FIGURE CAPTIONS 



Figure 1. The plot of X(r) versus r for three solutions of equation (40) with the boundary 
conditions that at r = 1, X = 2 and X r = —1,0,1. The function —X{r) is plotted in this 
figure as well. Note that the solution for X r (r = 1) = —1 approaches the special solution 
(41), also depicted here, for r — > 0. The behavior of the solution for X r (r = 1) = 1 near 
ro — 0.22 is in accordance with equation (43). 

Figure 2. The functions p(r)/Zo and <f>(r) given by equations (59) and (60), respectively, 
with 0o = tt/2 are plotted here as polar coordinates (i.e., abscissa = p cos0 and ordinate 
= p sin0). This represents the trajectory of a particle following the "radial" geodesic 
given by equations (58) — (61). The particle starts at r = from p = and moves 
counterclockwise until it returns to a singular axis (p = 0) at r = lo. The return trip takes 
only about ten percent of the total proper time Iq. 
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